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LDA + Negative U Solves a Puzzle of too Large Calculated Magnetic Moment in 

Iron-based Superconductor LaFeAsOi^Fa. 

H. Nakamura/'B N. Hayashi, 1 - 2 ^ N. Nakai, 1 ' 2 -@ and M. Machida 1 - 2 -| 

1 CCSE, Japan Atomic Energy Agency, 6-9-3 Higashi-Ueno, Taito-ku Tokyo 110-0015, Japan 
2 CREST (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan 
(Dated: June 30, 2008) 

A puzzle in the iron-based superconductor LaFeAsOi-^F^ is that the magnetic moment obtained 
by first-principle electronic structure calculations is unexpectedly much larger than the experimen- 
tally observed one. For example, the calculated value is ~ 2.0/iB in the mother compound, while 
it is ~ 0.3/xb in experiments. We find that the puzzle is solved within the framework LDA + U 
by expanding the U value into a slightly negative range. We show U dependence of the obtained 
magnetic moment in both the undoped x — 0.0 and doped x = O.f 25. These results reveal that the 
magnetic moment is drastically reduced when entering to the slightly negative range of U. Moreover, 
the negative U well explains other measurement data, e.g., lattice constants and electronic DOS at 
the Fermi level. We discuss possible origins of the negative U in these compounds. 

PACS numbers: 74.25.Ha, 74.25.Jb, 74.70.-b 



Very recently, iron based novel superconductors have 
attracted much attention, since its discovery may open 
another pathway toward the room temperature super- 
conductor. Consequently, the superconducting critical 
temperature has exceeded 50K within just a few months 
since the start of the high-T c race At the same time, 
several theoretical mechanisms have been proposed Q • 

Here, we summarize the experimental results accumu- 
lated since the discovery. The typical undoped mother 
compound LaFeAsO shows an antiferromagnetic spin 
density wave (SDW) ordering after the lattice structural 
transition from the tetragonal to the orthorhombic struc- 
ture at ~ 150K when decreasing the temperature d, 0). 
The SDW disappears with doping carriers via chemical 
substitution of a part of O by F, and the superconducting 
phase instead emerges from x ~ 0.05 in LaFeAsOi^^Fj;. 
Thus, the superconductivity has a close relationship to 
the magnetism and the lattice distortion. In addition, 
the optical absorption, photoemission and NMR stud- 
ies reported that a large pseudo-gap 19 meV ~ 100 meV 
opens even above the superconducting transition T c [f| . 
These imply an anomalous DOS suppression at the Fermi 
level by development of an order or other reasons. How- 
ever, the clear origin still remains unsolved. 

Besides the above measurement activities, the new ma- 
terial discovery has stimulated first-principle electronic 
structure calculations @, H, 0, E3] • Initially, inconsistent 
results on the magnetic structure were suggested. But, 
there is now a good agreement between the calculations 
and the experiments in respect of the magnetic and the 
related lattice structures. The antiferromagnetic SDW is 
confirmed to be the lowest energy state, in which the or- 
thorhombic structure is stabilized by optimizing the lat- 
tice structure under the SDW magnetic order [7[. How- 
ever, it is noted that there remains a big puzzle in the 
magnetic moment on an Fe atom in the SDW state. The 
moments calculated with the local-density approxima- 



tion (LDA) and the generalized gradient approximation 
(GGA) are ~ 1.0 Mb (LDA) to - 2.0 /x B (GGA) @, which 
is much larger than experimental data 0.2 ~ 0.35/ib ob- 
tained by the powder neutron scattering [3[, Mossbauer 
effect 11 1, and muon spin relaxation [l2j]. Normally, the 
moment over-estimation in LDA(GGA) calculation is al- 
most rare. This clearly implies that the calculations lack 
an unknown effect peculiar to the compounds or requires 
a particular correction. For example, we immediately no- 
tice that spin fluctuations always suppress the moment 
and two-dimensionality due to layered materials some- 
what enhances the fluctuation. However, S ~ 2 is too 
large to be reduced to about ~ 1/5 times only by the 
fluctuations. Thus, this puzzle may be crucial for eluci- 
dating the superconducting mechanism. In this paper, 
we report how the puzzle is solved within the framework 
of LDA + U. We believe that the result gives a hint on 
the superconducting mechanism. 

The Coulomb repulsion effect inside the local Fe d- 
orbit is listed as a missing effect. The effect is approxi- 
mately treated within "LDA + U" , U(= U e s) of which 
is theoretically decomposed into [/ e ff — U — J , where 
U is the on-site Coulomb repulsion (Hubbard U) and J 
is the atomic-orbital intra-exchange energy (Hund's pa- 
rameter) in the "simplified" framework [13J. We vary the 
single parameter U c s — U ~ J from a large positive value 
to a slightly negative one and also calculate electronic 
structures in another "separable" framework for U and 
J 1J| in order to identify which parameter is more essen- 
tial. We find that a slight negative U e g well explains the 
experimental results. 

The on-site Hubbard U employed in the electronic 
structure calculations is usually a positive value. The 
positiveness promotes the localized character of d- 
electrons and enhances the magnetic moment in the cases 
of magnetically ordered compounds. Since the parame- 
ter U is treated within the mean-field level in the stan- 
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FIG. 1: The crystal and magnetic structure in the lowest- 
energy state of the mother compound LaFeAsO calculated in 
the standard LDA (U e a = 0). The arrows on FeAs layers 
denote the spin direction. 



dard LDA+ U, more sophisticated methods coupled with 
DMFT and QMC have been also suggested [HGj] ■ How- 
ever, we point out that one cannot expect a drastic mo- 
ment reduction irrespective of the use of such advanced 
methods as long as the positive large U is taken into ac- 
count. In contrast, we extend the parameter U e g range 
to a negative one. Although the negative case is not pop- 
ular, it occurs in the following two cases. The Hubbard U 
itself is negative, and the intra-exchange J is effectively 
larger than the Hubbard U. The former case has been 
suggested by many authors in various theoretical mod- 
els as discussed later. Besides such explicit reasons, a 
slightly negative Hubbard U simply may compensate the 
energy calculation errors in the LDA calculations. How- 
ever, there has been just a few works taking the nega- 
tive on-site Hubbard U into account within the LDA+ U 

On the other hand, 



framework to our knowledge [171 . 
the latter case |20| may be also possible in the present 
compounds since J is suggested to be large while U 
is reported to be effectively smaller than our naive ex- 
pectation [2l], [23] • These facts require that one flexibly 
chooses U e ff in the present compounds at least. 

In this study, we employ the VASP code[23| in which 
the projector augmented- wave (PAW) method 24, 25| is 
implemented, and the lattice and the electronic struc- 
tures are optimized by choosing GGA + U. The sta- 
bilities of some different magnetic structures are exam- 
ined by comparing their total energy. The typical lowest- 
energy state at U = for the mother compound is shown 
in Fig. [TJ whose magnetic structure is equivalent with 
Ref. For the mother compound (x — 0), we actually 
adopt the La4Fe4As4C>4 cell that have space group sym- 
metry Ibam, in order to realize SDW antiferromagnetic 
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FIG. 2: The C c ff dependence for the undoped compound 
(solid curve) LaFeAsO and doped compound (dashed curve) 
LaFeAsO i-x Fa; (x = 0.125) on (a) the total energy difference 
between the SDW (see Fig. [TJ and the non-magnetic state, 
and (b) the SDW magnetic moment on a Fe-site, in which 
the experimental data (dotted curve) taken from Ref.0] is 
plotted. 



state as in Fig. [TJ In the case of the doped compound 
(x = 0.125), we use the super cell as shown in Fig. [TJ 
where one of O atoms is replaced by F. In our calcu- 
lation, the energy cut-off for the plane-wave set is set 
to 500 eV, the convergence condition for electronic self- 
consistent loop is less than 10~ 5 eV in the total energy 
difference, and the structure relaxation loop is repeated 
until all forces on ions are smaller than 0.02 eV/A. The 
grids for k points are taken as 6 x 6 x 6 and 4x4x2 for 
the undoped and doped compounds, respectively. 

Let us show numerical calculation results. Firstly, 
we focus on the mother compound LaFeAsO. Figure 
(2J]a) shows U c s dependence of the energy difference be- 
tween the non-magnetic and the stable SDW states, and 
Fig. [2](b) displays U c g dependence of the magnetic mo- 
ment in the SDW state. At C4ff = where several cal- 
culations were previously made, the moment is ~ 2.0 
Hb and the energy difference is ~ 90 meV per LaFeAsO 
cell. These values are consistent with the previous results 
0, 0, B E3] , in which the lattice structure is optimized 
and the orthorhombic one is observed to be stable. For 
U c s > 0, the effect of U e s further develops the moment 
value. As expected, the positive U c s stresses the local- 
ized character of d-orbital electrons, and overemphasizes 
the inconsistency with the experiments. On the other 
hand, the moment decreases with decreasing U c g and co- 
incides with the experimental results in a negative U e g 
region. In addition, we note that the energy difference 
becomes very small around U c g = —1, as in Fig. HJa). 
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TABLE I: The comparison between U = -1.0(J = 0.0) and 
U — 0.0(J = —1.0) on the magnetic moment and lattice con- 
stants for the mother compound in the separable LDA+ U. 



Structure [/^b] a [A] b [A] c [A] 



U — J 



U=-l J = 



U = J = 1 



Experiment" 



SDW 
NM 
SDW 
NM 
SDW 
NM 
SDW 



0.49 5.69166 5.67323 8.66216 

5.68683 5.68703 8.64482 

0.50 5.69166 5.67323 8.66216 

5.68693 5.68693 8.64482 

0.77 5.70995 5.65622 8.68032 

5.68693 5.68693 8.64482 

0.35 5.71043 5.68262 8.71964 



"Experimental data are taken from Refs. 0, Q] 



This indicates that SDW state becomes not so strongly 
stable as that at U e g = 0. 

Figure [3ja) displays U c g dependence of lattice con- 
stants, a, b, and c for the mother compound. The exper- 
imental values Q are shown by horizontal dashed lines. 
In this paper, b is the direction along which the antiferro- 
magnetic ordering grows. In the positive U e g range, the 
increase of U e g brings separations from the experimen- 
tal results for a and 6, while they show good agreements 
with experimental ones in the slightly negative range. 

From the results as seen in Fig. [5]and[3£a), it is found 
that the slightly negative U c g well explains the exper- 
imental results. Here, we check which parameter (U 
or J) is more essential in reproducing the experimen- 
tal results by using the separable framework [14|. Ta- 
ble U is a comparison between U = — 1.0(J = 0) and 
J = 1.0(U = 0) for the magnetic moment and the lat- 
tice constants. These results indicate that U = —1.0 is 
slightly closer to the experiments but it is not a conclu- 
sive difference. Moreover, the lattice constants obtained 
by the optimization do not differ significantly in both 
cases. We point out that a clear-cut determination is 
impossible here within the present framework. 

Next, let us turn to the doped compound. The cal- 
culations are made on the supercell as shown in Fig. 1, 
where one of O atoms is replaced by F and x = 0.125 
in LaFeAsOi-^Fa;. At the doping value, the compound 
shows the superconducting ground state. The dashed 
lines of Fig. [2jb) show U c g dependence for the magnetic 
moment value in the doped compound. The data cal- 
culated at U c g = still shows the stability of the SDW 
state, which is equivalent with the previous literatures 
0, @|. Then, the energy difference ( ~ 90 meV) is too 
large to replace the SDW state by the non-magnetic or- 
dered one even in the doped one. Also, the increase of 
U e g in the positive range enhances the stability of the 
SDW state, which is inconsistent with the experimental 
results. On the other hand, the negative U c g diminishes 
the moment resulting in the disappearance of the SDW, 
and the system recovers to the tetragonal lattice struc- 
ture by the optimization. The behavior of the lattice 
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FIG. 3: The I7 c ff dependence for (a) the undoped compound 
LaFeAsO and (b) the doped compound LaFeAsOi-zF^ (x = 
0.125) on the lattice constants a, b, c, and their ratio (a/b) 
in which the experimental data taken from Ref. are also 
plotted. The solid curve and dashed curve stand for the SDW 
and non-magnetic states, respectively. 



constants is also similar to the undoped case as seen in 
Fig. EJb). At the negative range, the calculated values 
except for c show good agreements with experimental 
ones. 

Figure [4] shows U c g dependent features of the elec- 
tronic DOS around the Fermi level in the case of the 
mother compound. Three cases, namely, U c g positive, 
zero, and negative ones are given in Figs.[4la), 4(b), and 
4(c), respectively. At U c g = 0, the gap like DOS sup- 
pression whose width is over ~ 300 meV opens, and the 
positive U c g simply expands the width, which results in 
much overestimation compared to the pseudo-gap value 
reported by recent experiments On the other hand, 

the slightly negative U e g gives the gap value much closer 
to the experimental ones (~ 30 to 100 meV) as shown in 
Fig. He). 

Finally, let us discuss the mechanism of the nega- 
tive U c g. Firstly, we give a summary of the studies on 
the negative Hubbard U . Many authors have suggested 
their ideas on the negative Hubbard [/, which are sum- 
marized in terms of theoretical target materials as fol- 
lows, i) several inorganic compounds forming the charge 
density wave, ii) amorphous semiconductors like chalco- 
genide glasses, iii) conducting polymers, and iv) heavy- 
fcrmion systems. These were intensively debated in old 
literatures related to the local electron pairing |26j |. On 
the other hand, more advanced studies inspired by the 
discovery of cuprate High-T c superconductors have been 
published on several metal oxides. A famous idea orig- 
inates from the charge disproportionation on the cation 
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FIG. 4: The electronic density of states DOS for (a) U e g = 2, 
(b) U e s = 0, and (c) U e « = — 1 in the undoped compound. 
E — corresponds to the Fermi energy. 



sites, which has been mainly discussed in doped super- 
conductors based on BaBiC>3, in which CDW occurs close 
to the emergence of its superconductivity [27j |. In ad- 
dition, the overscreening effects on the Coulomb inter- 
action coupled with the strong correlation (due to the 
low carieer density) have been proposed in the context of 
high-T c superconductivity mechanism [jjjj . If the over- 
screening, i.e., the attractive interaction effectively works 
between two electrons, then unusual softening of elec- 
tronically coupled longitudinal optical (LO) phonon com- 
pared to transverse optical (TO) one has been predicted. 
In fact, there are several reports which confirm the LO- 
TO frequency inversion due to the drastic LO softening 
in high-T c superconductor as well as other various metal 
oxides 29] . Thus, the negative Hubbard U is now not 
rare. But, since the effect is counterintuitive, more exper- 
imental and theoretical tasks are required. On the other 
hand, the idea that J is relatively effective in fixing U e /f 
may be very simple and acceptable. In this iron based 
superconductors, the Hubbard U may be rather small 
because Fe related five bands are entangled within their 
wide band-width. In fact, there are some experimental 
reports which support it (2lL l22| . 

In conclusion, we calculated electronic structure 
of iron based typical superconducting compound 
LaFeAsO^Fi-a; using the framework LDA+U with ex- 
panding the range of U(= U e g = U — J) from the pos- 
itive to the slightly negative range. Consequently, we 
found that the calculated magnetic moment, the lattice 



constants, and the pseudo-gap feature shows good agree- 
ments with the experimental results in the negative U e g 
range. We discussed some negative origin of U e g in the 
superconductor. 
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